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Investigation of the Cubic Boron Nitride Nucleation under the High 
Pressure and the High Temperature 
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In this paper we have theoretically found the activation energy(420.38fc J/mol) for the transformation from 
hBN to cBN in the microscopic viewpoint. We have introduced an analytical formula representing the 
dependence of nucleus formation time on the activation energy, synthesis pressure and temperature. We have 
theoretically determined the boundary line of cBN nucleus formation region in the P — T phase diagram on 
the basis of the diffusion mechanism of cBN nucleus formation. We have found that the cBN crystal nucleus 
formation time is less than 300?ns by comparing of the theory and experiment. 
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I. INTRODUCTION 

Cubic boron nitride (cBN) is usually produced in a 
high pressure-high temperature process by using mostly 
the hexagonal boron nitride(hBN) as the starting mate¬ 
rial. 

After the first discover of the cubic boron nitride syn¬ 
thesis method, so far, many studies on the cubic boron ni¬ 
tride synthesis was progressed but the mechanism of the 
synthesis reaction is still under discussion. Most exper¬ 
imental results described so fajlUE 

suffer from one clear 

disadvantage. 

Since all experiments for cBN synthesis were carried 
out in a closed high pressure chamber, they were not ob¬ 
served in situ. Only the final products after quench and 
pressure release were analyzed and only these results can 
be taken into account to model the transition mecha¬ 
nism. But the final process of pressure and temperature 
release may have an important influence. If the prod¬ 
ucts can only be analyzed after experiment, metastable 
phases may occur or vanish or other phases could un¬ 
dergo phase transition thus leading to incorrect conclu¬ 
sions about the appropriate mechanism. These disad¬ 
vantages are avoided if the observations are carried out 
under in situ conditions. In the papei^by using the syn¬ 
chrotron radiation for energy dispersive X-ray diffraction 
the cBN synthesis process was observed in situ. It is 
found that cBN grows directly from an amorphous ma¬ 
trix and not from different BN polymorphs as previously 
proposed. This is the evidence that the cBN formation 
process is carried out by the diffusion mechanism. 

In this paper, on the basis of these experimental facts, 
the diffusion mechanism of cBN synthesis were theoreti¬ 
cally studied and the cBN synthesis region in the P — T 
phase diagram was determined. 

II. THE CALCULATION OF THE ACTIVATION ENERGY 

The individual atom of hBN amorphous phase moves 
separately to form sp 3 bond with cBN crystal particle 


by diffusion, which would be overcome the energy bar- 
rier(activation energy). 

P+{n ) and P-(n) that are called the probability for 
jumping of an atom from hBN particle to cBN particle 
containing the n atoms and leaving of an atom from cBN 
particle to hBN particle are represented as follow^] 

y-y / \ ~ ( u 1 dAu\ 

— ■— j, (1) 

where h and uj are the number and frequency of atoms 
in the surface of cBN particle respectively, A G{n) is the 
change of the Gibbs free energy by forming cBN crystal 
nucleus containing n atoms, k is Boltzman constant, U is 
the height of the energy barrier that would be overcame 
for jumping of an atom from hBN to cBN. 

From now, let us estimate the activation energy, height 
of energy barrier that would be overcome for an atom 
to break away from sp 2 state of hBN to infinitely far 
dwelling site. 

The activation energy U was calculated using expres¬ 
sion (|2j) 

U = E 0 -E 1 -(E b + E n )/ 2, (2) 

where Eq and E\ are the energies of hBN perfect crystal 
and crystal having an vacancy occurred by the diffusion 
of atom, Eb and En are the spinpolarized energies for 
the isolated boron and nitrogen atoms. For calculating 
the energies of the isolated boron and nitrogen atoms, 
bcc crystal lattice of which B atoms are placed in the 
apices and a N atom is placed in the middle of a cube 
with 2.5 nm of edge length is introduced and local density 
approximation (LDA) is used. The calculation of energy 
of crystal having an vacancy is performed by using 16 
atoms and relaxation process. The precision of energy is 
O.OOOleV. 

The activation energy, calculated by using abinit ap¬ 
plication program according to the above method, is 

U = 420.38k J /mol. (3) 



2 


III. THE KINETIC THEORY OF CBN CRYSTAL 
NUCLEUS FORMATION ON THE BASIS OF THE 
DIFFUSION MECHANISM 

Now let us consider the kinetic theory of cBN crystal 
nucleus formation on the basis of the diffusion mecha¬ 
nism. Let us denote the cBN crystal nucleus number 
consisting of n atoms in t moment by f(n,t). Then the 
total number of cBN crystal nucleus is as following 


N{P,T,t) = Y,fM- ( 4 ) 

n =0 

The dependence of the f(n, t) on the time must satisfy 
the following equation 

df{n,t) 

— -Q t — = fi n ~ 1 ) t)P+{ n ~ 1 ) + 

f(n + l,t)P-(n + 1) - f(n,t)[P+(n) + P-(n)\, (5) 

where P±{n) is the probability for jumping of an atom 
from hBN to cBN and leaving of an atom from cBN to 
hBN given by formula ([l]) respectively. If n 1, equa¬ 
tion^) is approximated as following 

d -^> = A{[P_(n)-P + (n)]/(n,t)} + 

^ • ^2 {[-P+W+ P-(n)]/(n,t)}. (6) 


For quantitatively discussing on the kinetic theory of 
the cBN crystal nucleus formation, we must solve equa- 
tion(12) under a given initial and boundary conditions. 
Since the crystal nucleus does not exist in initial time, we 
can have the initial condition as following f(n, 0) = 0. In 
the moment of the crystal nucleus formation with radius 
of r, the Gibbs free energy change is as following 


AG = — r 3 n 0 (/x c - Hh) + 47rr 2 (cr c -ah), (13) 

where /i c and fih are the chemical potential of the cBN 
and hBN while a c and ah are the surface free energy 
density of cBN and hBN. no is the atomic number density 
of the cBN crystal. Since the crystal nucleus with small 
size must generate with big probability at a condition 
r —> 0 by the thermodynamic fluctuation theory, we can 
take the boundary condition as follows 


lim 

r —>-0 


/M) 

fe(r) 


—► o? 


(14) 


where f e (r) is the distribution function of the equilibrium 
state and is as follows 


fe(r) 



(15) 


It is the Fokker-Planck equation. However, in generally 
since U(n) can be approximated as following 

p -( n ) - p +{n) = nue-'Kr • (?) 

P+(n) + P-(n ) = 2fioje~^ r (8) 

P{n) = nwe-KT. (9) 

Using (7), (8)and (9), equation (6) is became as following 


where AG is the Gibbs free energy change given by for- 
mula(13). 

Then the distribution function f e (r) increase infinitely 
in out a critical size but the distribution function f(r,t) 
finite so that at case of the condition r —> oo we can take 
the boundary condition as following 


lim 

r—> oo 


/(u t) 

fe(r) 


0. 


(16) 


df{n,t) = 
dt 

J_d_ 
kT dn 


d' 2 
dn 2 


[P(n)f{n,t)\ + 


P(n)f(n,t) 


dAG 


dn 


In the case n 1, since 


dP{n) 

dn 


<C P(n) 


( 10 ) 


( 11 ) 


Let us solve equation (12) on cBN crystal nucleus of the 
critical size neighborhood under these initial and bound¬ 
ary conditions. In the case that we consider a problem in 
the critical size neighborhood of cBN crystal nucleus, we 
can expand AG'(n) in the neighborhood of critical size 
Uk by Taylor series and restrict to a second term. Since 
most of the non-steady state are near to the stationary 
state, generally, for the rough estimate of the non-steady 
time we can exchange approximately as following 


consequently, equation(lO) is became as follows 


df{n,t) = d_ 
dt dn 

_L_S_ ' 

kT dn 


p (n ) 


df(n,t ) 


P{n)f{n,t) 


dn 
dAG 


dn 


d_ 

dn 





d/oo 

dn 


d_ 

dn 


C Pf )-+ 


d_ 

dn 


(Pfoo) , 


(17) 

(18) 


( 12 ) 


where is the distribution function of the stationary 
state that reach through the non-steady process. This 
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stationary state distribution function is given by solving 
the stationary state Fokkcr-Planck equation as following 


d 

dn 



1 d 
kT dn 



dAG\ 
dn J 


0. (19) 


If we transform the equation (12) in consideration of the 
above facts then the equation (12) is transformed to fol¬ 
lowing simple form 


dj_ 

dt 


P kd 2 AG 
kT ^ dn 2 


(f-foo). 

n=rik 


( 20 ) 


Solving the equation (20) under initial and boundary 
conditions, we can obtain the following distribution func¬ 
tion of cBN crystal nucleus with critical radius 


/ = fco (l - e ~ t/T ) , (21) 

where /oo is the stationary state distribution function 
and t is as following 

18lnl0akTeVr . 

TraQUJnl(AV) 2 (P — a — bT ) 2 ’ 

where ao is radius of the volume space occupied by an 
atom, AV = V), — V c is the atomic volume difference be¬ 
tween hBN and cBN, a = —lA6GPa,b = 0.0031 GPa/K 
are constants determining the phase equilibrium line be¬ 
tween cBN and hBN and ui is the oscillation frequency 
of an atom in the cBN crystal surface, k is Boltzmann 
constant, P is the pressure ,T is the temperature, U is 
activation energy, a is the surface energy density be¬ 
tween cBN crystal nucleus and hBN phase and it can 
place a = cr c — <Th- As shown from formula (21), since 
the cBN crystal nucleus distribution function is arrived 
to the stationary state after r time,r is called non-steady 
time and we can call this time as crystal nucleus forma¬ 
tion time. 


IV. THE DETERMINATION OF THE CBN CRYSTAL 
NUCLEUS FORMATION REGION 

Representing the formula (22) in terms of the pressure, 
the following formula is obtained. 


P = a + bT 


3(2??rlO) 1 / 2 (cr/cT) 1 / 2 e2fcT 

(7rwT) 1 / 2 aQ|noA'i;| 


(23) 


Using the formula (231, we can find the cBN formation 


region in the P—T diagram. As known from formula (23), 
the P—T region in which cBN crystal nucleus is formed is 
determined by the activation energy and nucleus forma¬ 
tion time. For the activation energy we have used expres¬ 
sion (2). And the other material constants were equal as 



T/K 


Figure 1. cBN nucleation region by the direct transition 
experiment 11 ^ 1 - theory, 2 - experiment, 3 - hBN-cBN phase 
equilibrium line 


following^. ao = 1.12 • 10 10 m, no = 1.7 • 10 29 /m 3 , V c = 
5.9 ■ 10 _3 °m 3 , V h = 9.05 ■ 10- 3O m 3 , <r c = 4.72 J/m 2 , a h = 
3.40 J/m 2 ,w ~ 2.0 • 10 14 /s. 

Now we can obtain the boundary line of cBN nucleus 
formation region in P — T diagram by using formula (22) 
according to change of the crystal nucleus formation time 
r . However, if t is 300ms, theoretical line(curve 1 of 
Figure 1) coincides with experimental line (curve 2 of 
Figure 1). Therefore we can estimate cBN crystal nucleus 
formation time as less than 300ms. 


V. CONCLUSION 

In this paper the activation energy for the transforma¬ 
tion from hBN to cBN was theoretically found from mi¬ 
croscopic viewpoint and estimated to 420.38fcJ = mol. 
We have introduced an analytical formula representing 
the dependence of cBN nucleus formation time on the 
activation energy, synthesis pressure and temperature. 
Using this formula, the boundary line of cBN nucleus 
formation region in the P—T phase diagram theoreti¬ 
cally was determined. We have found that the crystal 
nucleus formation time of cBN is less than 300ms by 
comparing of our theoretical calculation and experiment. 
In this paper the temperature effect was not considered. 
If the temperature effect is considered, then the activa¬ 
tion energy will be decreased, so that the crystal nucleus 
formation time of cBN will be decreased. 
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